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Abstract 



A beam energy scan has been carried out at the Relativistic Heavy Ion Collider at Brookhaven National Laboratory to 
"^vQ -search for the onset of deconfinement and a possible critical point where the transition from a Quark Gluon Plasma to a 
y—i ^'hadronic phase changes from a rapid cross-over to a first order phase transition. Anisotropy in the azimuthal distribution 
of produced particles such as the second and third harmonics v 2 and V3 are expected to be sensitive to the existence of 
a Quark Gluon Plasma phase and the Equation of State of the system. For this reason, they are of great experimental 
interests. In this Letter we report on calculations of v% and V3 from the AMPT model in the Default(Def.) and String 



from 7.7 GeV to 2.76 TeV. 



^-s Mclting(SM) mode to provide a reference for the energy dependence of v 2 and V3 for 

'We expect that in the case that collisions cease to produce QGP at lower colliding energies, data will deviate from the 
AMPT String Melting calculations and come in better agreement with the Default calculations. 
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1. Introduction 

The motivation for colliding heavy ions at facilities like 
the Relativistic Heavy Ion Collider (RHIC) at Brookhaven 
National Laboratory and the Large Hadron Collider (LHC) 
at CERN is to form a state of matter called the Quark 
Gluon Plasma (QGP) Q. Each of these collisions creates a 
region so hot and dense that quarks and gluons become the 
relevant degrees of freedom instead of hadrons 0. Study- 
ing the conversion of coordinate space anisotropics into 
momentum space anisotropics gives insight into the nature 
of the matter created in these collisions [3j- For decades, 
elliptic flow (v 2 — (cos 2(0 — ^rp))) has been studied to 
probe the conversion of the elliptic shape of the initial over- 
lap zone into azimuthal anisotropy in momentum space 0] 
over a broad range of colliding beam energies. Measuring 
the strength of that conversion as a function of beam en- 
ergy to search for evidence of the onset of deconfinement 
or a softening of the equation-of-state is one of the goals 
of the RHIC Beam Energy Scan program. In 2007 Mishra 
et. al. [B| proposed the analysis of v\ for all values of n 
and argued that density inhomogeneities in the initial state 
would lead to non-zero v„ values for higher harmonics in- 
cluding Vs. Although they proposed that v n vs. n could be 
used to search for superhorizon fluctuations, it was later 
noted that higher harmonics of v n would be washed out 
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by viscous effects and that the shape of %, vs. n would 
provide a valuable tool for studying 77/s 6]. It was also 
subsequently pointed out that information on was to a 
large extent contained within already existing two-particle 
correlations data [7j , and that v n and v n fluctuations would 
provide a natural explanation for the novel features seen 
in those correlations, such as the ridge like [8j and mach- 
cone like [9( structures. That the ridge could be related 
to flux-tube like structures in the initial state was already 
argued by Voloshin in 2006 [loj . Calculations carried out 
within the NEXSPHERIO model showed that in a hydro- 
dynamic model fluctuations in the initial conditions lead 
to a near-side ridge correlation and a mach-cone like struc- 
ture on the away-side 11] . In 2010, Alver and Roland used 
a generalization of participant eccentricity (£n,part) to arbi- 
trary values of n as in Ref. [H and showed that within the 
AMPT model, the final momentum space anisotropy for V3 
was proportional to the initial £3 jPa rt [HI- This explained 
the previous observation that the AMPT model produced 
correlations similar to those seen in the data (albeit with 
smaller amplitudes) [3]. Later studies showed that with 
changes to the input parameters, AMPT could quantita- 
tively describe the centrality dependence of v 2 and ^3 at 
200 GeV and 2.76 TeV Q3- 

In this paper we use the AMPT model to study the 
beam energy dependence of v 2 and V3 . We study collisions 
ranging from y/s NN — 7.7 GeV to 2.76 TeV. We compare 
AMPT in the string melting and the default setting. 
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Figure 1: The charged particle multiplicity density scaled by 
iVpart/2 in the AMPT model for String Melting and Default modes. 
The red line shows the parameterization of experimental data pre- 
sented in Ref 1191 . 



2. AMPT Settings and Comparisons to Data 

The AMPT model provides two modes: Default and 
String Melting AMPT in default mode is essentially 
a string and minijets model (without a QGP phase) where 
initial strings and minijets are produced with the HIJING 
event generator (r7l |. The interactions of the minijet par- 
tons are then calculated using a parton cascade (ZPC) be- 
fore the strings and partons are recombined and the strings 
are fragmented via Lund string fragmentation. ART (A 
Relativistic Transport model for hadrons) is used to de- 
scribe how the produced hadrons will interact. In the 
String Melting mode, the strings produced from HIJING 
are decomposed into partons which are fed into the par- 
ton cascade along with the minijet partons. The partonic 
matter is then turned into hadrons through coalescence 
and the hadronic interactions are subsequently modeled 
using ART. So while the Default mode describes the colli- 
sion evolution in terms of strings and minijets followed by 
string fragmentation, the String Melting mode includes a 
fully partonic QGP phase that hadronizes through quark 
coalescence. The model therefore provides a convenient 
way to investigate expectations for a variety of observables 
with and without a QGP phase. 

Several parameters need to be specified in the model 
including parameters a and b for Lund string fragmen- 
tation, the QCD coupling constant a s (which the model 
treats as a constant), and the screening mass for gluons 
in the QGP phase \x. A recent study found that a good 
description of the multiplicity density, v 2 and V3 could be 
achieved with the parameter set: a=0.5, b=0.9 (GeV -2 ), 
a s =0.33 and ^=3.2 (fm _1 ) [l5|. In this study, we found 
that we can acheive a good desciption of the multiplicity 
density at all energies from ^Ts^= 7.7 GeV to 2.76 TeV by 
using parameter set: a=2.2, b=0.5 (GeV -2 ), a s —0A7 and 



(M=1.8 (fm _1 ) and turning off initial and final state radia- 
tion in HIJING. In this case, the initial cutoff for minijets 
Po does not need to be adjusted with y/s in order to match 



the LHC multiplicity densities [18f. We leave po and all 
other parameters fixed for all energies. Figure [T] shows 
the charged particle multiplicity density scaled by N part /2 
for 0-5% central Au+Au or Pb+Pb collisions from AMPT 
String Melting and Default vs The line shows the 

parameterization of the experimental data from Ref. (l9j . 
Both the SM and Default calculations are in good agree- 
ment with the experimental data throughout the entire 
energy range. 
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Figure 2: Elliptic flow data from AMPT and experiments at 
v /s^-=62.4 GeV, 200 GeV [STAR], and 2.76 TeV [ALICE]. For the 
String Melting calculation we show V2 calculated relative to the par- 
ticipant plane «2{PP} defined by the positions of the nucleons and 
using the two particle cumulant t>2{2} = (cos2(cAi — <f>j)). Experi- 
mental results are shown for the two-particle V2{2} and four-particle 
i>2{4} cumulants. 
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In Figure [2l the AMPT model results are compared 
to experimental data at ^/s^ =62.4 GeV, 200 GeV, and 
2.76 TeV. For the SM calculations, we show (i) vi rela- 
tive to the participant plane calculated from the initial 
conditions of AMPT and (ii) the two-particle cumulant 
results f2{2}. The difference between those results rep- 
resents the amount of non-flow correlations present in the 
model. All the model calculations have a similar centrality 
dependence but the Default results are well below the SM 
results. The data generally agree well with the SM calcula- 
tions. In the case that the v 2 fluctuations are dominated by 
eccentricity fluctuations, i>2{4} should be equal to Vi with 
respect to the reaction plane [20| . The fact that the exper- 
imental «2{4} results are slightly below the model results 
calculated with respect to the participant plane does not 
therefore signify a discrepancy between data and model. 
We consider the agreement between the model and the 
data to be satisfactory. 
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Figure 3: Top: v 2 (pr) calculated with AMPT SM relative to the 
reaction plane (^{RP}) for beam energies from 7.7 GeV to 2.76 TeV. 
Bottom: The t;2{RP} data at different energies are shown scaled by 
the results at 200 GeV. 

STAR has shown that for p T < 1 GeV, v 2 {4}(p T ) in- 
creases with y/s NN , for pt > 1 GeV «2{4}(j»t) is roughly 
independent of collision energy in the range 7.7 GeV to 
2.76 TeV [21]. It is surprising for a measurement that is 
supposed to be sensitive to viscosity and collective effects 
in the expansion to not depend on y?s NN over such a wide 
range of energies where the initial conditions and proper- 
ties of the fireball should be changing quite significantly. 
Given this surprising experimental result, it is interesting 
to see if the same trend is reproduced in the AMPT model. 
In Figure HI we show V2{pt) calculated with respect to the 
reaction plane for collisions with center of mass energies 



ranging from 7.7 GeV to 2.76 TeV. Although the statistics 
in our study were not sufficient to calculate V2{4:}(pt), it 
has been shown that as long as v 2 fluctuations are dom- 
inated by eccentricity fluctuations and those eccentricity 
fluctuations are Gaussian distributed along the x and y 
axis, then «2{4} is equivalent to «2{RP} [20j]. We there- 
fore check to see if v 2 {RP} (pt) is independent of y/s NN 
for pt > 1 GeV in the AMPT model. We find that the 
variation of t> 2 {RP} is not large in AMPT throughout the 
energy range studied. For px < 1 GeV, 1>2{RP} varies 
by about 5% from 7.7 GeV up to 200 GeV. Going from 
200 GeV to 2.76 TeV, v 2 {RP} increases by 20%, inde- 
pendent of pr. In the RHIC range, the AMPT w 2 {RP} 
results for pt > 1 GeV are actually increasing as the en- 
ergy is decreased with w 2 {RP} at Pt = 1-5 GeV for 7.7 
GeV being 20% larger than for 200 GeV. This likely re- 
flects the softening of the spectrum which allows flow ef- 
fects that push low momentum particles to higher momen- 
tum, to have a larger influence at intermediate Pt- The 
same trends hold when studying v 2 {PP}(pt) (not shown). 
Although there are differences between the trends seen in 
AMPT and in the data, one can conclude that even in the 
AMPT model, the changes in v 2 {RP}(p T ) or v 2 {PP}(p T ) 
when increasing ^/s^ from 7.7 GeV to 2.76 TeV are not 
large. In this case, it is not necessarily surprising that 
the data also does not change drastically. Since based on 
the AMPT model, we would not expect a large variation 
of V2{BP}{pt) with y/s NN , as long as one assumes that 
a string melting or QGP phase exists throughout the en- 
ergy range under study, the fact that the data seem to 
change very little no longer appears to be so difficult to 
understand. 

3. The Third Harmonic 
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Figure 4: The slope of (113} vs. £3 as a function of the square root 
of the number of participants for four different colliding energies. 

Having shown that our parameter selection provides a 
good description of the charged particle multiplicity den- 
sities and the elliptic flow, we now turn to investigate V3 
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Figure 5: t>3{2} and «3{PP} from AMPT SM and Default calcualtions for y'sjvjv = from 7.7 GeV to 2.76 TeV. Experimental results are 
shown at 2.76 TeV fjl. 



and its energy dependence. We first study the relation- 
ship of v 3 to the third harmonic participant eccentricity. 
In Ref. the AMPT model is used to show that v 2 
and v 3 have a linear dependendence on e-i and £3. At 
Quark Matter 20f f , STAR showed that v 3 je 3 scales with 



1 1 v/iVpart [22| • Here we check to see if this phenomenolog- 
ical observation is also reproduced in the AMPT model. 

In Figured] we investigate the dependence of the slope 
of (V3) vs. £3 on N part - The figure shows d(v 3 )/de 3 vs. 
yHV^Tt for =7.7 GeV, 39 GeV, 200 GeV, and 

2.76 TeV. We find that for all the energies investigated 
(including those not shown in the figure), d(v 3 }/de 3 in- 
creases linearly with y/N pal f The AMPT model there- 
fore correctly describes the phenomenological observation 
made by STAR. This also indicates that according to the 
string melting version of AMPT, even at energies as low 



as i/s]v7V=7.7 GeV, v 3 reflects the fluctuations in the ini- 
tial geometry of the collisions and that the centrality de- 
pendence will remain similar at all energies although the 
magnitude will change. At the lowest energies investigated 
here, the contributions from jets and minijets should be 
negligible so they will not contribute significantly to the 
centrality dependence of v 3 . The experimental observation 
of a similar centrality dependence for v 3 at 7.7 and 200 
GeV 22[, therefore strongly contradicts assertions that v 3 
is dominated by jet-like correlations [23[ . 

In Figure El AMPT SM and Default calculations of 
v 3 {2} and v 3 {PP} are shown for 7.7, 11.5, 19.6, 27, 39, 
62.4, 200 GeV and 2.76 TeV. While v 3 {PP} reflects the 
true correlation of particles with the initial participant 
plane, «3{2} includes non-flow effects. The difference be- 
tween v 3 {2} and i> 3 {PP} is large at 200 and 39 GeV while 
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at 7.7 GeV v 3 {2} and u 3 {PP} are equivalent. This indi- 
cates that indeed, according to AMPT SM, non-flow docs 
not make an appreciable contribution to «3{2} at 7.7 GeV. 
We compare the model results to ALICE data at 2.76 TeV 
and find that w 3 {PP} for AMPT SM matches the AL- 
ICE data on v 3 {2}. The v 3 {2} AMPT SM results over 
predict the ALICE data and the v 3 {PP} AMPT Default 
results underpredict the ALICE data. The t>3{2} Default 
results also underpredict the ALICE data for N part > 100. 
The correspondence of u 3 {PP} from AMPT SM with v 3 {2} 
from ALICE data means that either non-flow is overpre- 
dicted in AMPT or v 3 is underpredicted. The 200 GeV 
data is also in good agreement with preliminary STAR 
data [22| (not shown) in the same centrality range. In more 
peripheral collisions, the STAR data in Ref. [22] tends to 
increase as seen with the AMPT «3{2} results. This sug- 
gests that while «3{2} measurements for -/V part > 100 are 
dominated by the correlation of particles with the partic- 
ipant plane, in more peripheral collisions i>3{2} begins to 
reflect correlations related to mini-jet structure similar to 
that in p + p collisions. 
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Figure 6: The 



dependence of v? j {2} (SM), v%{PP} (SM), 



and w|{PP} (Default) for two different centrality intervals. 

In figure [5] we show the AMPT results for the variation 
of v%{2} and v|{PP} with ^/J^ from 7.7 GeV to 2.76 TeV 
for two centrality intervals. The results on u|{PP} using 
the default setting for AMPT are very small and well be- 
low the preliminary data presented by STAR. The u|{PP} 
SM results decrease rather smoothly with decreasing en- 
ergy but still have an appreciable value down to 7.7 GeV. 



The calculations for u|{2} SM have the same value as the 
i;f{PP} SM at 7.7 and 11.5 GeV. This again indicates that 
within this model, nonflow from minijets has a negligible 
impact on two-particle correlations at the lowest energies 
measured in the RHIC beam energy scan. Above those 
energies, the difference between i>|{2} and v|{PP} grows 
substantially. It will be interesting to see if the experi- 
mental data on v 3 follows the same trend as AMPT SM 
all the way down to 7.7 GeV where non-flow from mini- 
jets can be neglected. It will be most interesting to see 
if data eventually drops down to the values predicted by 
the AMPT Default model. Estimates of the Bjorken en- 
ergy density [25| compared to the Lattice QCD estimates 
for the critical energy density [26| sug gest that this may 
not happen until below 7.7 GeV [27|. The calculations 
presented in this paper provide a base-line with which to 
compare future experimental data. 



4. Viscosity Estimates 

In Ref. [28[ , the authors estimate the shear viscosity to 
entropy ratio rj s / s based on kinetic theory using 



Vb/s 



3tt 



1 



40a 



^9 + f^ 2 



(1) 



where /i is the screening mass of a gluon in the QGP and 
a s is the QCD coupling constant. Both are input param- 
eters for the AMPT model. In Ref. [HI , it is shown that 
parameter set B (/1 = 3.2 fm _1 and a s — 0.33) provides 
a good description of the V2 and v 3 data at y/s NN = 200 
GeV and 2.76 TeV. These parameters and equation Q] then 
yield the rj s /s(T) as shown in figure [7] with the label Set 
B. For the estimated initial temperatures (378 MeV) at 
top RHIC energy, r] s /s=0.38, far above the ADS/CFT or 
quantum lower bound of 1 /Air [29], |3(| and also above most 
estimates from hydrodynamic models. 

In this study, we have used parameter set A (^ = 1.8 
fm _1 and a s — 0.47) corresponding to a larger partonic 
scattering cross section of 10 mb rather than 1.5 mb for 
set B. We find a good description of the data with set A 
by turning of initial and final state radiation in HIJING. 
This reduces the initial entropy production and multiplic- 
ity but in a way that matches the multiplicity at all the 
energies studied without varying any other parameters. 
The smaller initial multiplicity compensates for the larger 
cross sections so that the data are still well described. The 
r] s /s(T) estimated from this parameter set is labeled Set 
A in figure For set A, r] s /s at 378 MeV is 0.088 which 
is very close to the ADS/CFT conjectured lower bound. 
We conclude therefore that the AMPT model can give a 
good description of the v% and v 3 data with a wide range of 
r] s /s values and that it's crucial to understand the initial 
entropy production in order to extract the correct value of 
rj s /s in the QGP phase. 



5 




200 250 300 350 400 450 500 
T [MeV] 

Figure 7: The temperature dependence of r) s /s according to eq. [l] 
for two parameter sets. 



5. Summary 

In this paper, we have presented AMPT SM and De- 
fault calculations for v 2 and U3. We found that we can 
describe RHIC and LHC data on multiplicity, v 2 and V3 
by turning off initial and final state radiation in HIJING 
(reducing the initial entropy) but keeping relatively large 
cross-sections in the QGP phase. These settings give a 
good description of event multiplicities from y/s NN =7.7 
GeV up to 2.76 TeV. Whereas a previous study found 
a good description of data using a much smaller cross- 
section implying a much larger value for rj s /s, our studies 
with a larger cross-section implies a ratio of viscosity to 
entropy much closer to the ADS/CFT conjectured lower 
bound. We have also studied how V2{pt) changes from 7.7 
GeV to 2.76 TeV. We find that within this model, v 2 (pt) 
changes very little across the whole energy range studied, 
consistent with what is observed in data. We also find 
that AMPT reproduces the experimental observation that 
^3/63 oc J V par t . Our studies of the centrality and beam 
energy dependence of v 2 and V3 with SM and Default set- 
tings provide a comparitive base-line for studies of v 2 and 
V3 in the RHIC beam energy scan. 
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